Role of d-wave pairing in ^415 superconductors. 
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We argue that the recent Raman spectroscopy observation of a sharp mode in s-wave supercon- 
ducting VaSi is due to a competing d-wave pairing state. We present microscopic arguments for the 
origin of this d-wave order. We further argue that the d-wave order explains the observed shrinking 
^— ^ , of the vortex core structure at anomalously low magnetic fields and the large anisotropy observed 

in the upper critical field. 

Competing pairing channels have played a prominent role in theories of superfluid 3 He, cuprate, and pnictide 
' superconductors. In p-wave 3 He, an excited pair state, stemming from a competing /-wave channel, has been observed 
through interferometry in an acoustic cavity [H . The observation of this /-wave pairing channel leads to predictions of 
rich physics near impurities, surfaces, and vortices. In the cuprates, a competing s-wave or d xy order has been invoked 
to interpret the properties of surface Andreev bound states [2|. More recently in the superconducting pnictides, many 
theories find that a sign changing s-wave state and a competing d-wave state have close transition temperatures [3|, 
^ . leading to the prediction of a pair collective mode observable through Raman spectroscopy and a superconducting 
5— i ' state that breaks time reversal symmetry 

In this work we argue that A 15 superconductors provide an ideal candidate for competing superconducting order. 



ON 



We show it is natural to expect an attractive d-wave pairing channel in addition to the attractive s-wave pairing 
channel in Alb materials. Specifically, the d-wave pairing belongs to the cubic E g representation. This representation 
has two degenerate states with the same symmetry properties as k% — k g and {2k 2 z — k% — ky)/\f?>. The primary 
Q ■ experimental motivation for exploringthis state is Raman spectroscopy. Raman active superconducting gap modes 
have been seen in NbsSn and in VaSi [5]. Modes have been found in three symmetries: Ai g , E g , and F2 g . The modes 
in the Eg channel are significantly sharper and occur at lower energies than those in the A\ g and F2 g channels. We 
argue that the E g modes are due to an attractive d-wave pairing channel that arises from the unique structure of 
O . cubic Alb materials. This d-wave order does not play an important role in the ground state thermodynamics but, as 

we show below, it is important for understanding the superconducting state of A15 materials in magnetic fields. 
— i • Many A15s are type II superconductors with transition temperature (T c ) on the order of 20 K 6]. Their cubic 
structure (space group Pm3n) contains the formulae unit A3B. The B atoms form a cubic lattice while the A atoms 
form three sets of orthogonal chains that run through the faces of the cube [2j. These chains play an important role 
in understanding these materials. For example, interactions on these chains drive the Martensitic cubic to tetragonal 



■ transition that appears in some of these materials 0, Q . Early theories assumed the electronic states on the chains to 
be one-dimensional 1(J 11|- However, it was later found that this was an oversimplification fl2j [l3l|. Nevertheless, the 



band structure shows strong anisotropics that originate from electrons residing on the chains [141 ]. In accord with this, 
y—i I we assume that the superconducting condensation energy is dominated by electrons on the three orthogonal chains. 

This naturally leads to three order parameters \& = (ip x ,ipy,ip z ) corresponding to superconductivity on each of the 

three chains (we call this the chain basis). In the following we exploit this chain basis to show how a d-wave collective 
• »h , mode arises that can account for the Raman spectroscopy data. We then turn to the development of a microscopic 

model that provides an estimate for the pairing interactions in the s and d-wave channels. Finally, we consider the 
jj^ 1 appearance of the d-wave order when magnetic fields are applied. 

Ground state and d-wave collective mode: To address the origin of collective modes, we begin with a Ginzburg Landau 

(GL) theory that is based solely on symmetry arguments. We use the chain basis to write this theory, in which case 

the GL free energy density is 

/ = aj2 i^i 2 + + \ D ^\ 2 + K2 J2 \ D ^\ 2 + k *J2 WiXA^)* + « 4 J2 

+ P n J2^ + ^ + cc)m 2 ) (1) 
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FIG. 1. Position of Vanadium atoms in the unit cell and the corresponding hopping matrix elements included in the Hamiltonian. 

where Di — —idi — 2eAi/c, A is the vector potential, and c.c. represents complex conjugate. Depending upon the 
coefficients e and Eq. Q] allow for four possible ground states: ^ s oc (1,1,1), 'Fd oc (1, e j27r / 3 , e l47r / 3 ), ^d,i oc 
(1,-1,0), or ^d,2 oc (—1,-1,2). The state ^> s has s-wave symmetry (A% g ) while the states V^d, ^d.i, ^^,2 all are d- 
wave, belonging to the E g representation. The state ^ s will be stable if e < while one of the states <Fd, ^di, or 1 $>d2 
will be stable if e > 0. Physically, it is expected that gap on each chain is maximized to maximize the condensation 
energy. This implies that the states ^d.i and ^,2 are unlikely to be realized. Of the remaining two states, the 
most likely is ty s as the preponderance of evidence suggests an s-wave ground state in ^415 superconductors. In the 
following it is assumed that ^ s is the ground state. However, it should be noted that it is possible that ^d is the 
ground state of some Alb superconductors. 

We now turn to possible collective modes of Since our order parameter has six complex degrees of freedom, 
there can exist six collective modes: the amplitude and phase modes for each of the chain gaps. The amplitude modes 
are likely to be strongly damped since the mode gaps will typically lie near the chain gap values. Furthermore, the 
phase mode in which all three chain order parameters oscillate in phase will be pushed to the Plasma frequency. This 
leaves two phase modes as candidates for collective modes. These modes correspond to relative phase oscillations 
between the chain order parameters and resemble Leggett modes of multiband superconductors [15j. To examine 
these two phase modes in more detail, we follow an approach developed in Ref. |16| in the context of unconventional 
superconductors. We represent the order parameter as \& = ipo(e l91 , e 1 ® 2 , e t6:s )/y/3. We use the Josephson relation 
hd8 a /dt — -~n a together with the quasi-conservation laws hdfi a /dt = (d/j, a / 'dN a )(dF / d6 a ) where fi a is the chemical 
potential, F is the free energy, and N a the particle number corresponding to quasi-particles on chain a. This yields 
a two-fold degenerate relative phase mode with d-wave E g symmetry that has a frequency given by 



[3|e|Vo 2 - (4/3 2 + fii + 3& + 2/3 6 )V 4 ] 



(2) 



where N(0) = dN a /dfj, a , is the density of states for each of the chains. 

In the limit e = = At = /3s = Pa = 0, this gives uj = 0. In this limit, the chains are uncoupled and the theory 
has an accidental U(l) x U(l) x U(l) symmetry associated with the independent variations of the phase of each chain. 
The microscopic theory presented later reveals that 4/3 2 + (3 4 + 3/?5 + 2/3 6 < 0, so that u is well defined. If cj < 2A , 
where A is the chain gap value, these modes will be weakly damped and provide a natural explanation for the Raman 
spectroscopy results of Ref. @ . 

Microscopic theory: Now we turn to a microscopic model to examine how an attractive ci-wave pairing channel can 
appear in A15 materials and also to determine the phenomenological coefficients in the GL theory. LDA calculations 
show that near the Fermi energy, the 3z 2 — r 2 orbitals on the A atoms that form chains along the z-direction and 
symmetry equivalent orbitals for the other two chains directions are the most important [l7^ . We therefore use a 
tight binding Hamiltonian that includes a nearest neighbor (nn) intrachain hopping (ti) as well as a next nearest 
neighbor (nnn) interchain hopping (i 2 ) between these orbitals (as shown in Fig. [TJ. This is similar to the approach 
used earlier by Bhatt [IH . To identify the interactions that lead to superconductivity, we note that there is compelling 
evidence for strong electron phonon (ep) coupling in A15 materials. In particular, phonon anomalies due to strong ep 
coupling of the chain atom displacements to electrons have been observed [l9| that are in agreement with predictions 
from LDA calculations 17]. Additionally, tunneling data on superconducting NbaSn reveal features of the phonon 
spectrum [2 Qj , showing that this coupling is responsible for superconductivity. The phonons that couple most strongly 



FIG. 2. The self-consistent d x 2_ y 2 like solution for the E g symmetry gap function shown on the Fermi surface (for fe/ti = 0.22). 
The gap values are represented by the different colors show in the box on the right. The (k x , k y ,k z ) are in units of the inverse 
lattice spacing. 



to electrons are intra-chain phonons jljjj ]. We therefore use the following BCS Hamiltonian as a minimal model for 
A15 superconductors 

H = t ^Y Y c li<r c »i° + * 2 Y Y c lia c ^3<r-VoY c in c ti c ivlc i ^, (3) 

VO <ij> nn VV' (7 <ij>nnn W 

where c y i a destroys an electron with spin er, with the Wannier function centered on the i th Vanadium atom that 
transforms as the v orbital [v corresponds to one of (3z 2 — r 2 , 3a; 2 — r 2 , 3y 2 — r 2 )]. The tight binding Hamiltonian is 



given in more detail in the supplemental notes 27J . The pairing interaction Vo corresponds to an effective interaction 
that would arise from ep coupling on the chains alone. The ep coupling gives rise to a local interaction that is retarded 
in time by an amount At ~ 1/uid, where Ud refers to the Debye frequency. This retardation is captured within BCS 
theory by restricting the energies of the electrons to lie within a cutoff ujd of the Fermi energy. 

To proceed, we diagonalize the kinetic energy term and this leads to three Fermi surface sheets and a momentum 
dependent pairing interaction V a ,p{k, k) (a, (3 correspond to the different bands). Treating this interaction within 
a mean-field theory and defining a fe-dependent gap on each of sheets of the Fermi surface, we get the following 
linearized gap equation 



A„(fc) = ln (Ij^r) ^(^/ 3 (fc,fc')A / ,(fc')^(fc')>F.5 (4) 

V Bey ~ 

where a,j3 correspond to the different bands, ep(k ) is the band energy, to = nuin, Np{k') — fe'[ep(k')]\ i s the 
density of states evaluated at the Fermi energy in the gap equation and the angular bracket implies average over the 
Fermi surface. To solve this equation we do not assume a fc-dependence for the gap but rather discretise the Fermi 
surface and explicitly find the eigenfunctions that solve this gap equation. We find that there are only two solutions, 
an s-wave A\ g solution and a d-wave E g solution. It is pos sible to include Coulomb interactions analytically in the 



limit t^fti = by using an Anderson- Morel type theory [211 ] ■ In this limit, we find that on-site Coulomb interactions 
have the same effect on both the s and d-wave solutions. 

In Fig. 2, we show the self-consistent solution for the E g gap degree of freedom corresponding to a d x i_ y i symmetry. 
This solution is qualitatively different than the /c 2 — k g solution expected for a spherical Fermi surface. In fact, it 
is more appropriately described as having three nearly constant gap amplitudes (denoted by A, 0, —A) which can 
be interpreted as the three chain gap values. In Fig. 3, we show the effective interactions for the A\ g and E g 
pairing channels as a function of til^x while keeping the Fermi surface volume constant (in the weak-coupling limit 
T c = u>d exp (^-) where A is the effective interaction found from solving the gap equation). 

Magnetic Fields: A15 superconductors exhibit anomalous properties in magnetic fields. Notably, these cubic materials 
have a large anisotropy in the upper critical field (H C 2) [22| and they also exhibit a hexagonal to square vortex lattice 
transition at anomalously low magnetic fields [23| . This vortex transition is also accompanied by a sharp reduction in 



the coherence length 2J, |25| . Here we examine whether the appearance of d-wave order can provide an explanation 
for these observations. 
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FIG. 3. Ratio of eigenvalues A^/As vs tijt\. 



First we begin with the anisotropy in the upper critical field. In particular, H c2 in V3S1 shows an anisotropy that 
is estimated to be 20 % at the zero temperature [22}. H c2 is highest for the field along the (100) direction, lower for 
fields along the (110) direction, and lowest for the field along the (111) direction .22]. As shown in the supplemental 



notes 



271 ] , the GL theory of Eq. [T] above predicts that the same anisotropy in H C 2 as observed experimentally (as 
described in the previous sentence) . Furthermore, this conclusion is independent of the values of the coefficients m . 

Now we turn to the hexagonal to square vortex lattice transition that occurs at low fields [23[. The free energy 
of Eq. Q] contains too many undetermined coefficients to carry this out. We therefore use the weak-coupling th eory 
described above to find the coefficients in the GL free energy. The details are given in the supplemental notes 27 1. 
Of particular relevance, we find that the coefficients /32,/?5, k&, and k% are relatively small and we neglect them in the 
following. The GL theory can then be used to construct a London theory valid at low fields. We take ipi = ipoe 1 ^, 
where \ipo\ takes on the bulk uniform solution value. This yields the following London free energy density 



h = fo + (e|^o| 2 + /34|Vo| 4 )[cos(^ - 2 ) + 



cos 



- 03 ) + COs(0 3 - 01 ) - 3] 



Ki|V»o| a 53(Vi& - 2eA if + K 2|^o| 2 X}(V0j - 2eA) 



(5) 



where /o is the free energy in the homogeneous state. Let us define 0i = 6 + 771 — /)a/v3, 02 = — r\\ — r/2/y/3, and 
03 = 9 + 2?72/-\/3. We expand /l to quadratic order in 771 and 772 (this is valid when Ki/e « a 2 , where a is the 
spacing between vortices). We then solve for 9, r\\ and 772 to find the following non-local London free energy 



h 



8- 



[B 2 + X 2 (W x B) 



4«i 



-{d x d y B z 



(6) 



3|e| Aq 3k 2 + K\ 

where Ao is the penetration depth. Eq. [5] has the same form as the non-local London theory used by Kogan [26[ 
to examine the hexagonal to square transition. The important point here is that the non-local correction (the third 
term in Eq. [6]) is due to appearance of d-wave order and it is possible for this term to become large at relatively 
low fields, when Ki/e«a 2 Ri -B/$ (^0 is the superconducting flux quantum). This would push the hexagonal to 
square transition to low fields, in agreement with experiment. However when K\/t > a 2 , the non-local London theory 
no longer applies. To understand the physics in this regime requires numerical analysis. Some understanding can 
be found by setting e = and /34 1 V^o 1 2 = 0. In this case, the three order parameters ipi become decoupled and the 
vortex cores for each component are smaller than that of the pure s-wave theory. This is consistent with the observed 
shrinking of the vortex core [13, [III . While more analysis is required, the above shows that the low field hexagonal 
to square transition and shrinking of the vortex cores is consistent with the appearance of a d-wave component in 
magnetic fields. 

In conclusion we have argued that in A15 superconductors an attractive ci-wave pairing channel exists in addition 
to the s-wave pairing channel. This d-wave pairing channel accounts for the sharp Raman resonance seen in the 
s-wave superconducting state. Furthermore, this d-wave pairing is induced by magnetic fields and can account for 
the observed anisotropy in the upper critical field and the observation of a sharp reduction of the superconducting 
coherence length at a hexagonal to square vortex lattice transition that appears at anomalously low fields. 

This work was supported by NSF grant DMR-0906633. We thank Steve Kivelson, Sonny Rhim, Manfred Sigrist, 
Kazuo Ueda, and Michael Weinert for useful discussions. 
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TIGHT BINDING MODEL 

The unit cell for V3S1 contains six effective Vanadium atoms per unit cell. We use a tight binding model that 
included the (3z 2 — r 2 ) orbitals in the z-direction, and the corresponding orbitals along the x and y directions given 
by group rotation operations that satisfy cubic symmetry Density functional theory shows that these orbitals are the 
most important near the Fermi surafce The Vanadium atoms in a unit cell along the x, y, and z chain directions 
have position r m = (if/ 2 ± x/4, z/2 ± y/4, z/2 ± y/4) in units of lattice spacing with respect to the center of the 



unit cell and correspond to orbitals with symmetries v = (3x — r , 3y 



, 3z — r ) respectively. We label these 



1 through 6 as shown in Fig. 1. We include nearest (intra-chain) and next-nearest neighbor (inter-chain) hopping 



J_V p -ik-{R t +r„ 



><? , where dL { annihilates an electron on the m 

rnK m i % 



th 



matrix elements (see Fig. 1). Letting <? m i — 
Vanadium atom that has the symmetry of the v orbital in the unit cell given by the Ri, the tight binding Hamiltonian 
in reciprocal space is expressed in matrix form H = ^2j. Mk^h with 



and 4< 



k 



( c ifc>' 



•2fc' u 3fc' u 4fe 



' C 4k > C ! 



1 

5ft' C 6kJ 



( -A* 

e x (k) 

Sl4(*) 
-15 (*) 

V 4 6 (k) 



e x (k) 

4 S (k) 

e* 4 (fc) 



ei3(fe) eu(fe) eis(fe) eie(fc) \ 

£23(fc) £24(fc) £25(fc) £26(fc) 



-A* 
e»(fc) 



^46 (*) 



£35 (fc) £36 (fc) 

£45 (fe) £46 (k) 

-H e y (k) 

e y (k) -fx 



(1) 



In the above 6x6 matrix \i is the chemical potential that is adjusted to yield 



the correct qualitative form for the known Fermi surface sheets as evaluated from the full DFT calculations 2] . Here 
the terms e a (k) = 2t\ cos(fc a /2) where a = {x,y, z} and the interchain terms are given by 

£ 24 (fe) = 2t 2 e- l[k - /i) e-' 1 ^/^ cos(k y /2) 
e 23 (k) = 2t 2 e-^!^e l ^l^ cos(k y /2) 
e u (k) = 2t 2 e l ^ /4) e~ l ^/ 4 '> cos(k y /2) 
e 13 (fc) = 2t 2 e 2( ^ /4) e' i( ^ /4) cos(fc a /2) 
e 26 (k) = 2A 2 e- 2(fc * /4) e-^/ 4 ) cos(fc z /2) 
e 25 (k) = 2t 2 e- l ^ /4) e l ^y^^ cos(fc z /2) 
e 16 (fc) = 2t 2 e l{k * /A] e- l{k y /A] cos(fc z /2) 
eis(fc) = 2t 2 e l ^ /4) e^ k y^ cos(fc,/2) 
e 46 (fc) = 2t 2 e- t(k ^^e-^ k y /A) cos(k x /2) 
e 45 (fc) = 2t 2 e- l{k */^e l{ - k y'^ cos(fc x /2) 
e 36 (fe) = 2A 2 e 4 ( fe ^ 4 ) e - l( ' £! ' /4) cos(k x /2) 

e 35 (k) = 2t 2 e^/^e i(k y^ cos(k x /2). (2) 



DERIVATION OF GINZBURG LANDAU THEORY 



In addition to the tight binding Hamiltonian above, a superconducting instability is present in the system through 
the following quartic interaction term, 



2 




Hbcs 



-Vq 



E 



(c 7 



m,i,\ 



771,1,4--' m,*,4, m,*,'I" 



(3) 



For simplicity the BCS Hamiltonian has been approximated as an onsite interaction which is a reasonable approxi- 
mation owing to the intra-chain phonons being the most strongly coupled phononsQ. 

Fourier transforming and diagonalizing the band part of the Hamiltonian gives a momentum dependance to the 
quartic interaction. In terms of the eigenstates of the band Hamiltonian, the full Hamiltonian can be expressed in the 
momentum space, 

H= J2 ^ a (k)ai(k)a a (k) +Y / V a ^(k,k')al(k + q/2,t)aU-k + q/2,i)a p (-k' + q/2,i)a p (k' + q/2A) (4) 

a=[l->6] a,/3 

where a, /? represent the eigenstates of the diagonalized band Hamiltonian. If c = P^a, where U = P* is the transpose 
of the unitary eigenvector matrix whose components are Uij then the momentum dependance of the quartic interaction 
term is 

v a Ak,k') = -VoY,K a ( k M a (- k nA-k')u t Ak'). (5) 

i 

The gap function is A Q ^(fc, q) — J2k' p Va,p(k, k ')( c k+q/2 p t c -k+q/2 p j.) wnere a i P represent band indices and 
inhomogeneity of the order parameter is included through the center if mass momentum index q (this will allow the 
GL coefficients ft, to be determined). The superconducting gap equation becomes 

A a {k,q)=T > — 6) 

o-u ( w « + e M fc ) + l A <3( fc )l) 

n,P,K 

where ui n = (2n + l)nT + iq ■ Vu'Cp(k'). 

Within a weak coupling approximation the gap equation close to the superconducting T c can be integrated to yield 
the Ginzburg Landau free energy. The general method for performing this operation is well knownQ. Here we assign 
the order parameters ip s to the A\ g representation solution and (771, 772) to the E g representation of cubic symmetry. 
The superconducting gap can be expressed in terms of these order parameters in the form, 

A(fc) = ^U(k)+ m f m (k)+ m f V2 (k) (7) 
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Here is the s-wave basis function and (f Vl , f V2 ) correspond to the basis functions for the two dimensional E ff 
irreducible representation of cubic symmetry. Expanding the gap equation to second order in q and third order in the 
order parameter yields an equation that can be integrated in terms of the order parameter to obtain the Ginzburg 
Landau free energy. Note that the basis functions f(k) are the self consistent numerical solutions of the linear gap 
equation. 

Therefore using the above procedure and expanding the superconducting gap equation in the weak coupling limit 
close to the superconducting transition temperature the free energy in the s-d basis is obtained, 

F = a a \i>\ 2 + «d|*?| 2 + <\DtJj\ 2 + k' 2 (\D Vi \ 2 + \D m \ 2 ) + K' 3 [(2{D z r) 2 )(D zm )* - (D x r l2 )(D xV2 y - (D y r) 2 )(D y r) 2 )*) 
-{2{D z r ll ){D z r ll y - (D xm )(D xm )* - (D ym )(D ym y) + V3((D xrh )(D x r] 2 y - (D ym )(D yV2 y + c.c)] 

+n' i [^={2{D z ^){D z r l2 y - {D x ^)(D xm y - (D y ^){D ym y) + ((A^XA^i)* - (D y ^)(D ym )*) + c.c] 

+ &M 4 + &\v ■ v\ 2 + P'M'dml 2 + \m\ 2 ) + &W 2 (ri + vl) + c.c] 
+/3' 6 [r(m\m\ 2 - 2m\m\ 2 - vlvt) + c.c] (8) 

Note that this is the most general free energy for the given representation allowed by cubic symmetry up to the fourth 
order in the expansion parameters. The weak coupling fourth order coefficients are numerically evaluated from the 
following Fermi surface averages 

(2 r2\ 



P'2 = c(f 2 Jl)FS 

^3 = 2U 2 J 2 2 )fs 



^5 - ^Ufuf^FS 

P' 6 = c(flU)FS 

(9) 

here the constant c = 7VbC(3)/(87r 2 T c 2 ). Similarly the gradient term coefficients are given by 

K' x = c'mfy F ) FS (io) 

4 = c'/Hflv F ) FS (ii) 

^^c'/V3(f m Uv Fx )FS (12) 

K> , 4, = c'(Uf Vl V Fx )FS (13) 

(14) 

where v F is the fermi velocity. 

The corresponding free energy in the chain basis has been given in the text. The chain basis is formally given by 

i>x = tps/V^-ili/VQ + m/V^, 4> y = ?Ps/V3- rji/VE-r] 2 /y/2 and ip z = iJ) s /y/3-r)iJ~%. The free energy in the chain 

basis follows from Eq.[H]by the following relations: /3[ = f3' 2 = (fii + /3 3 )/3, /3 3 = /3' 5 /2 = /3' 6 /(2V2) = (2ft - /3 3 )/12, 
/3' 4 = (4/?! + /3 3 )/3, k[ = k' 2 = (3ki + 2k 2 )/5, k' 3 = -^JIk' a = -(3/10)(ki - n 2 ). The variation of the GL coefficient 
with the interchain chain hopping is shown in Fig. 2. Note that the largest coefficients are /3 3 , ki, and k 2 and it 
would be reasonable to simplify the theory including only these terms in a GL analysis (as we have done) . A similar 
simplification does not occur in the s-d basis. 



p2 r2\ 



UPPER CRITICAL FIELD ANISOTROPY 



The upper critical field may be evaluated by expressing the free energy in terms of the creation and annihilation 
operators which operate on the Landau levels. For a magnetic field along the z direction they are defined as, 

U+ = -j={D x +iD y ) (15) 
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FIG. 2. Coefficients of the GL theory vs interchain interaction. The volume enclosed by the Fermi surface has been kept 
constant with the variation in ta/ti- 
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V2 



(D X ~ iDy) 



2e 
Tic 



H 



(16) 



(17) 



where I 2 = jfff, with H being an external magnetic field applied along the z direction, and D a = —id a — 2eA a /c. 
Expressing the equivalent operators for the different directions of magnetic field we can then write the quadratic 
portion of the free energy for the magnetic field along the (1, 0, 0) and (1, 1, 0) direction in the matrix form as, 



/loo = ( i>* V* V2 ) 
/no = ( r Vi V*2 ) 



f KiIIo + a 



— K4II1 
K1II0 + a s 





^n 2 -« 4 ni 

K 2 n - k 3 U 2 + a d -V3K3II1 

—VSK3II1 K 2 n + K3II2 + a d 

W 2 

K 2 n - K3II2 + a d 

K 2 n + k 3 it 2 + a d 




(18) 



(19) 



(20) 



Here n = 1// 2 (1 + 2n+n_), III = l/(22 2 )(n + Ii\ + Hi), n 2 = l/(2/ 2 )(n - 3(n2_ + n 2 .)). Note that the (1,1,0) 
free energy has been written in the rotated basis with, 

x = z 
y' = {x- y)/V2 
z' = (x + y)/V2 

(21) 

Comparing the upper critical field along the (1,0,0) and (1,1,0) magnetic field directions, we find that the free 
energy for the (1,0,0) direction contains additional off diagonal terms which couple the order parameters of the 
two dimensional representation with <i-wave symmetry and another term which mixes the d-wave and s-wave order 
parameters. Since the presence of these off diagonal coupling terms will enhance the upper critical field in general, the 
(1, 0, 0) field direction will have a higher H c2 than the H c2 along (1, 1, 0). The role of s-d mixing in comparing the H c2 
for the (1,1,1) and the (1,1,0) magnetic field direction requires a more detailed approach. For \T — T c \ « T s — T4, we 



■5 



can use a perturbation theory which treats the Hamiltonian corresponding to the s-d mixing as a perturbation (with 
the s-wave state is the ground state). Defining the second order correction (E^) to the s-wave eigenvalue as 

l -^ = ±+E^ (22) 

we have 

E (2) = y f \<0\H>\ Xlm >\2 + | < 0\H'\x2m >\ 2 \ (23) 
\ E — E\ m E — E 2m J 

h[ w = h/(2i 2 ){i + 2n+n_ + 3(n 2 + n 2 )) (24) 
# in - VW+ + n 2 _) (25) 

(26) 

In the above expression for the perturbed Hamiltonian, the 2D representation of cubic symmetry has been expressed 
in the r\\ = (1,lj,uj 2 ), ry 2 = n{ basis, where w = e 27 ™/ 3 . Also the states \ Xm > are the eigenstates of the unperturbed 
Hamiltonian derived from the diagonalization of the 771,772 free energy. The solutions can be expressed in terms of 
Landau level wave function. For the field along the (1,1,0) direction the eigenstates are 

Xim = <Pm{eix,y/ei) (27) 
X2m = 4>m{t2X,y/e 2 ) (28) 

2 





- 2k 2 




+ «2 




-f 2k 2 


Kl 


- «2 



4 = \l^— - m 

where (f> m represent the m th Landau level solution. Similarly for the field along the (1,1,1) direction the eigenstates 
Xim, X2m are of the form, a m \m > +b m \m + 2 > where a m and b m are known functions. It is found upon solving this 
problem numerically that the second order correction to the energy will be given by, 

(2) _ 8(« 4 /«i) 2 r „ n 
E ^~~IHT S -T d ) (31) 

(2 ) _ 9.5(« 4 /«i) 2 . . 

^iio~ lHTs _ Td) W 

therefore the (1,1,0) direction to provides a higher upper critical field than the (1,1,1) magnetic field direction near 
T c . The GL theory given above thus yields the same anisotropy as seen experimentally, independent of the values of 
the coefficients Ki. 
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